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Abstract 

In this paper, numerical and solitonic solutions of Korteweg de Vries 
(KdV) and Korteweg de Vries-Burger's (KdVB) equations with initial 
and boundary conditions are calculated by sinc-collocation method. The 
basis of method is sine functions. First, discretizing time derivative of 
KdV and KdVB's equations using a classic finite difference formula and 
space derivatives by 9— weighted scheme between successive two time lev- 
els is applied, then Sine functions are used to solve these two equations. 
Mathematica programming is used to solve matrix representation of these 
equations. KdV equation describes behaviorof traveling waves which is a 
third order non-linear partial differential equation (PDE). Maximum ab- 
solute errors are given in Tables. The figures show approximate solutions 
of these two equations. Three conservation laws for KdV's equation are 
obtained. 

Keywords: Numerical method; KdV equation; KdV-Burger equation; Sine 
method; Collocation. 

1 Introduction 

First, In 1895, Korteveg and de Vries derived famous KdV equation that de- 
scribes weakly nonlinear shallow water waves and models one directional long 
water wave of small amplitude , propagating in a channel. It occurs in many 
field of physics such as in water waves, plasma and fiber optics. Another ex- 
ample for this equation is pulse wave propagation in blood vessels. After its 
discovery scientist found solution of this equation that is called soliton. The 
dynamic of solitary wave is modeled by this equation. 

KdVB equation was derived by Su and Gardner [1] for a wide class of nonlin- 
ear system in the weak non-linearity and long wavelength approximation. The 
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steady state solution of the KdVB equation has been shown to model [2] weak 
plasma shocks propagation perpendicularly to a magnetic field. When diffu- 
sion dominates dispersion the steady state solutions of the KdVB equation are 
monotonic shocks, and when dispersion dominates, the shocks are oscillatory. 
The KdVB equation has been obtained when including electron inertia effects 
in the description of weak nonlinear plasma waves [3] . The KdVB equation has 
also been used in a study of wave propagation through liquid field elastic tube 
[4] and for a description of shallow water waves on viscous fluid. 

Consider third order partial differential equations 

u t + 6uu x + u xxx = 0, x G ft = (a, b) C K, t > 0, (f ) 

as KdV equation. KdV's equation has the analytical solution as 

u(x, t) = 0.5sec/i 2 (0.5(a; - t)), (2) 

and consider 

Ut + euu x ~ vu xx + ptu xxx = 0, x e Q = (a, b) C M, t > 0, (3) 
as KdVB's equation that has analytical solution given as 

u{x > *) = w [1 + tanh( w {x + £ <)} - l sech2{ i {x + (4) 

that space variable is defined as 

Xi = a+(i-l)h, i = l,...,N, h= J^^l (5) 

as KdVB's equation. In recent year numerous methods are used for solving 
KdV and KdVB's equations. Wang Ju-Feng et al. obtained numerical solu- 
tion of the third-order nonlinear KdV equation using the elementfree Galerkin 
(EFG) method which is based on the moving least-squares approximation. A 
variational method is used to obtain discrete equations, and the essential bound- 
ary conditions are enforced by the penalty method [5] . The dynamics of solitary 
waves is modeled by the Korteweg de Vries (KdV) equation. Jamrud Aminuddin 
and Sehah, strated by discreetizing the KdV equation using the finite difference 
method. The discreet form of the KdV equation is put into a matrix form. The 
solution the of matrix is determined using the Gauss- Jordan method [6]. Jie 
Shen et al. have studied the eventual periodicity of solutions to the initial and 
boundary value problem for the KdV equation on a half-line and with periodic 
boundary data. They derived a representation formula for solutions to the lin- 
earized KdV equation and rigorously establish the eventual periodicity of these 
solutions [7]. Julio Duarte et al. employ the Wavelet- Petrov- Galerkin method to 
obtain the numerical solution of the equation Korterweg-de Vries (KdV) [8] . Ho- 
motopy Perturbation Method (HPM), Variational Iteration Method (VIM) and 
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Homotopy Analysis Method (HAM) for the semi analytical solution of Kortweg- 
de Vries (KdV) type equation are applied by Foad Saadi et al [9]. Also, KdVB's 
equation has been solved in recent years. M. T. Darvishi et al. have used a Nu- 
merical Solution of the KdV-Burgers Equation by Spectral Collocation Method 
and Darvishis Preconditionings [10]. In [11] Riccati equation expansion method 
is presented for constructing exact travelling wave solutions of nonlinear evolu- 
tion equations. The main idea of this method is to take full advantage of the 
more solutions of Riccati equation to construct exact travelling wave solutions 
of nonlinear evolution equations. More new exact travelling wave solutions are 
obtained for KdVB equation. On using variable transformations and proofs of 
theorems, the asymptotic behaviour and the proper analytical solution of the 
Korteweg-de Vries-Burgers equation have been found in [12]. Anna Gao et al. 
studied the problem of optimal control of the viscous KdVBs equation, they 
develop a technique to utilize the Colc-Hopf transformation to solve an optimal 
control problem for the viscous KdVBs equation [13]. In [14] exact travelling 
wave and solitary solutions for compound KdVBs equations are obtained by 
using an improved sine-cosine method and the Wu elimination method. Nu- 
merical solutions of the Korteweg-de Vriesequation using the periodic scattering 
transform /i- representation is studied in [15]. Operator Splitting Methods for 
Generalized KortewegDe Vries Equations has been discussed in [16]. 

The paper is organized into six sections. Section 2 outlines some of the main 
properties of sine function and sine method. In Section 3, the discretization 
of KdV and KdVB equations is discussed. Section 4 outlines stability analysis 
and section 5 introduces errors and conservation laws. Finally numerical results 
and the efficiency and accuracy of the proposed numerical scheme is shown by 
considering some numerical examples in Section 6. 



2 The Sine function 

In this section the basis of sine function is discussed [17]. The sine function is 
defined on the whole real line, — oo < x < oo, by 

Smc(x) = < 7TX 7-1 fQ\ 

w \ 1, x = 0. w 

For any h > 0, the translated sine functions with evenly spaced nodes are given 
as 

S(j,h)(z) = smc(^±), j = 0,±l,±2,-... (7) 

The sine functions are cardinal for the interpolating points Zk = kh in the 
sense that 

S(j,h)(kh)=6% = { l> k k = j j (8) 
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If / is defined on the real line, then for h > the series 



00 ' u 

C(f,h)(z)= J2 f(jh) S mc(^-), 



] = -co 



(9) 



is called the Whittaker cardinal expansion of / whenever this scries converges. 
They are based in the infinite strip D s in the complex plane 



D s = {w = u + iv : \v\ < d < — }. 



(10) 



Some derivatives of sine function will be used in reduction the equation to matrix 
form so [13], 



lj°> = [S(j,h)(x)]\ x=Xi 



1, j = i: 
0, 



(11) 



and 



_ 1 



0, 3 = i; 

(_l)(fc-3) 



7F1T 



1 

' 0, 



3 ' j — 

-2(-l)< fc - J '> 
(k-j) 



3 = 



(-1) 



(fe-i) 



And so on, for even coefficients, where r — 1, 2, . . . 



r(2r) 



[S(j,/i)(z) 



7r2r (-l) r 
/i 2r+l 



3 = v, 



(-1) (J !j V r ~ 1 r iV+l 2r! 2< /■-• „-\2i ■ /■ 



(12) 



(13) 



H6-7r 2 (fc-. 7 ) 2 L jVi, (14) 



(15) 



and for odd coefficients, where r — 1,2, ... 
d 2r 



r(2r+l) 



dx 2 ^ 1 



[s(j,ft)W] 







3 = «; 



. (-1)" 2 , y-lf iy (2r+l)l ot/j ,-\2J _• /,• 



(16) 
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3 Survey of the method 

Consider third order partial differential equations 

u t + 6uu x + u xxx = 0, x G ft = (a, b) C R, t > 0, (17) 

as KdV equation and 

Ut + euu x — vu xx + [iu xxx = 0, x G Q = (a, b) C K, t > 0, (18) 

as KdVB's equation, with the initial condition 

u(x,0) =f(x), xeTi, (19) 

and the boundary conditions 

u{a, t) = g a (t), u(b, t) = g b (t), t > 0, (20) 

where e, fi and f are constants. By discrediting time derivative of KdV's equa- 
tion using a classic finite difference formula and space derivatives by 9- weighted 
scheme we have 



u n+l _ u n 



+ 0(e(uu x T +1 +v(u xxx r +1 ) 



St 

+ (l - 6)(e(uu x ) n + ^(u xxx ) n = 0, (21) 

so using Taylor expansion for the term uu x and considering e = 6 and /j = lwe 
have 

u n+1 + Ste(6[u n u^ +1 + <w™ +1 ] + u n + x ) 

= u n + (6(29 - l)(uu x r + St(l - 6)u n xxx ) (22) 

for KdV equation and with same calculation for KdVB equation we have 



u n+1 + 8te(e[u n ul +1 u 1 iu n+1 ] - vu n x t x + pLu n + x ) 

= u n + 8tdu n u n x - 6t{l - 0)[euu n x - vu n xx + nu n xxx \ (23) 

Now we use approximate solution as 

N 

u(x,t n ) = u n (x) ~Y, u ? S j( x )- ( 24 ) 
i=i 

where 



Sj(x) = Sin( U h. )■ (25) 



,. r x — (j — l)h — a x 
~~h 

By substituting above approximate solution in Eq. (22) a matrix representation 
is obtained for KdV equation as 



Mu n+L = R (26) 
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where 

A d = ] ■ i = 2, ...,N- 1, j = 1, ...,iV, elsewhere] NxN , 

A b = [l\f : i = l,N,j = 1,..., N, elsewhere] NxN , 

B = -[I^p : i = 2,...,N- l,j = 1,...,N, elsewhere] NxNl 

C = [l\f ■ i = 2,...,N- l,j = 1,...,N, elsewhere] NxN , 

G = -[l\f : i = 2, . . . ,N - l,j = 1, . . . ,N, elsewhere] NxN , 

u™ = Bu n , D=p(u n ) p - 1 *u2*A d , E=(u n )*B, 

F n+1 ^[g a (t n + 1 ) 7 0,-.-,0,g b (t n+1 )} T , 

so with these definition for KdV equation we have 

M = [Aa + A b + 95t(6(E + D) + G)], 

R = [A d - 5t{6{26 - l)E - (1 - 8)G}]u n + F n+1 , 

and with same substitution for KdVB's equation 

M = [A d + A b + 95(s(E + D)+^G- vC)], 

R=[A d + 5t{e{26 - l)E - (1 - 9){jiG - vC)}]u n + F n+1 , 



4 Stability Analysis 

In this section stability analysis of approximate solution for linearized equation 
is discussed [18]. The error at nth time level is 



4.1 KdV's equation 

By considering the obtained matrix we have 

[H + 5t6K]e n+1 = [H- St(l - 6)K]e n , (27) 

where H = [A d + A^A- 1 and K = [6E + G]A~ 1 so, 

e n+i = pe n 

where P = [H + StOK]- 1 ^ - 6t(l - 6)K). This method is stable if || P || 2 < 1 
or p(P) < 1 which is spectral radius of the matrix P. The stability is assured 
if all the eigenvalues of the matrix [H + StOK]^ 1 ^ — St(l — 9)K] satisfy the 
following condition 

X H - 5t(l - 6)\ K 



x h + stex K 



< 1 (28) 
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where Xh and Xk are eigenvalues of the matrices H and K respectively. When 
6 = 0.5, the inequality (28) becomes 



X H -0.55tX K 



X H + 0.56tX K 



< 1. 



(29) 



In the case of complex eigenvalues Xh = o-h + «6^and Xk = o-k + ^fc, where 
Ofc, dk, bh and 6^ are any real numbers, the inequality (29) takes the following 
form, 

(ah — 0.5#tdfc) + i(bh — 0.5Stb k ) 



(a h - 0.5Sta k ) + i(b h - 0.55tb k ) ~ 

The inequality (30) is satishcd if a^a^ + b^b k > 0. For real eigenvalues, the 
inequality (29) holds true if either (Xh > and > 0) or (Xh < and A^ < 0). 
This shows that the scheme is unconditionally stable if a^afc + bhbk > 0, for 
complex eigenvalues and if either (Xh > and A^ > 0) or (Xh < and A^ < 0), 
for real eigenvalues. When 8 = 0, the inequality (29) becomes 



1 - 



5t\ k 



X h 



< 1, 



i.e., 



St < 



2A 



H 



and 



A 



H 



> 



Xk Xk 

Thus for 9 = 0, the scheme is conditionally stable. The stability of scheme for 
the other values of can be investigate in a similar manner. The stability of 
the scheme and conditioning of the component matrices H, K of the matrix P 
depend on the weight parameter and the minimum distance between any two 
collocation points h in the domain set [a, b]. 



4.2 KdVB's equation 

By considering the obtained matrix we have 

[H + 5t6K}e n+1 = [H- St(l - 6)K]e n , (31) 
where H = [A d + A b ]A^ and K = [eE + [iG - vC\A~ x so, 

e n+l = p e n 



where P = [H + StOK}- 1 ^ - St(l - 6)K}. This method is stable if || P || 2 < 1 
or p(P) < 1 which is spectral radius of the matrix P. The stability is assured 
if all the eigenvalues of the matrix [H + StOK] -1 ^ — St(l — 0)K] satisfy the 
following condition 

Xh - 6t(l - 9)X K 



A 



H 



stex 



K 



< 1 



(32) 
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where Xh and Xk are eigenvalues of the matrices H and K respectively. When 
6 = 0.5, the inequality (32) becomes 

X H -0.h8tX K 



X H + 0.55tX K 



< 1. 



(33) 



In the case of complex eigenvalues Xh = o-h + i&^and \ K = a k + ib k , where 
ah, dk, bh and b k are any real numbers, the inequality (33) takes the following 
form, 

(a h - 0.5Sta k ) + i(b h - 0.5Stb k ) 



(a h - 0.5Sta k ) + i(b h - 0.5Stb k ) " (34) 

The inequality (34) is satisfied if a^a k + bhb k > 0. For real eigenvalues, the 
inequality (33) holds true if either (Xh > and > 0) or (Xh < and < 0). 
This shows that the scheme is unconditionally stable if a^afc + bhb k > 0, for 
complex eigenvalues and if either (Xh > and A^ > 0) or (Xh < and Afe < 0), 
for real eigenvalues. When 8 = 0, the inequality (33) becomes 

5tX k 



1 - 



< 1, 



i.e., 



St < 



2X H 



and ^— > 
Xk 



Thus for 9 = 0, the scheme is conditionally stable. The stability of scheme for 
the other values of can be investigate in a similar manner. The stability of 
the scheme and conditioning of the component matrices H, K of the matrix P 
depend on the weight parameter and the minimum distance between any two 
collocation points h in the domain set [a, b]. 



5 Errors and Conservation Laws 



In this section, two error norms is defined that will be used for showing the 
accuracy of the method as follows 



L2 =|| u — u W2- 



Loo =11 u - u 



N 



h^2\u-Uj\ 2 , 
max\uj — Uj\, 1 < j < N 



3 3 1 

where u, u are exact and approximate solution respectively. 
KdV equation has three conservation laws as follows 

fb 



h 
h 
h 



u(x, t)dx, 



~dx, 



f u(x,tf 
\ [u(x,t) 2 - -u(x,t) 3 ]da 

J a 3 
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Where Ii, I2, I3 represents mass, momentum and energy that show conserva- 
tion properties of collocation method by, 

N 
N 

N 1 

/ 3 ^]T[K) 2 --K) 3 ], 

Sinc-collcation scheme satisfies the properties I\ and I2 in the conservative case. 



6 Numerical Solution 

In this section L 2 and are obtained and shown in Tables and approximate 
solution of KdV and KdVB's equations is shown in Figures. 

6.1 KdV equaiton 

By considering /j, = 1, e = 6, 6 = 0.5, in Table 1, two kinds of error is calculated 
for n = 100, a = -15, b = 15, St = 0.1, T = 0.1, . . . , 0.9. Three invariants of 
conservative laws is obtained. Table 2 indicates errors for St = 0.01 and Table 
3, shows error for deltat = 0.001. Table 4, errors are given for different time 
levels T = 0.001, . . . , 0.009. By Tables, Errors reduce when time step decreases. 
Figure 1, indicates solitonic solutions of KdV equation in different time level 
of T = 1,. . .,9 for n = 100, a = -10, b = 20, St = 0.01. Figure 2, indicates 
solitonic solutions of KdV equation in different time level of T = 1, ... ,9 for 
n = 100, a = -15, b = 15, St = 0.1. 



Time 




L 2 


h 


h 


h 


0.1 


4.55798x10" 


-b 


7.31340x10" 


-b 


2 


0.666667 


-0.0888889 


0.2 


9.28126x10" 


-5 


1.34140x10" 


-4 


2 


0.666667 


-0.0888853 


0.3 


1.23469x10" 


-4 


1.81748x10" 


-4 


2 


0.666667 


-0.0888766 


0.4 


1.34961x10" 


-4 


2.20735x10" 


-4 


2 


0.666667 


-0.0888665 


0.5 


1.39977x10" 


-4 


2.54695x10" 


-4 


2 


0.666667 


-0.0888567 


0.6 


1.58066x10" 


-4 


2.85701x10" 


-4 


2 


0.666667 


-0.0888476 


0.7 


1.77771x10" 


-4 


3.15070x10" 


-4 


2 


0.666667 


-0.0888394 


0.8 


1.98187x10" 


-4 


3.43619x10" 


-4 


2.00001 


0.666667 


-0.0888321 


0.9 


2.20096x10" 


-4 


3.71849x10" 


-4 


2.00001 


0.666667 


-0.0888257 



Tablel: Errors and invariants for 
n = 100, a = -15, b = 15, St = 0.1, T = 0.1, . . . , 0.9. 
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Time 


Loo 


L 2 


h 


h 


h 


0.1 


1.74342x10" 


b 


1.91251x10" 





2 


0.666667 


-0.0888889 


0.2 


2.31194x10" 


-6 


2.87693x10" 


6 


2 


0.666667 


-0.0888888 


0.3 


2.73488x10" 


-6 


3.64966x10" 


6 


2 


0.666667 


-0.0888887 


0.4 


3.04988x10" 


-6 


4.32025x10" 


6 


2 


0.666667 


-0.0888886 


0.5 


3.32879x10" 


-6 


4.95178x10" 


6 


2 


0.666667 


-0.0888885 


0.6 


3.49763x10" 


-6 


5.53733x10" 


6 


2 


0.666667 


-0.0888884 


0.7 


3.93719x10" 


-6 


6.16498x10" 


6 


2 


0.666667 


-0.0888883 


0.8 


1.98187x10" 


-6 


6.75605x10" 


6 


2.00001 


0.666667 


-0.0888883 


0.9 


4.69034x10" 


-6 


7.26574x10" 


6 


2.00001 


0.666667 


-0.0888882 



Table2: Errors and invariants for 
n = 100, a = -15, b = 15, St = 0.01, T = 0.1, . . . , 0.9. 



Time 


Loo 


L 2 


h 


h 


h 


0.1 


1.74199x10" 





1.75980x10" 


b 


2 


0.666667 


-0 


0888889 


0.2 


2.30437x10" 


6 


2.54006x10" 


-6 


2 


0.666667 


-0 


0888889 


0.3 


2.71876x10" 


6 


3.16315x10" 


-6 


2 


0.666667 


-0 


0888889 


0.4 


3.07014x10" 


6 


3.70810x10" 


-6 


2 


0.666667 


-0 


0888889 


0.5 


3.27326x10" 


6 


4.23733x10" 


-6 


2 


0.666667 


-0 


0888889 


0.6 


3.69103x10" 


6 


4.73099x10" 


-6 


2 


0.666667 


-0 


0888889 


0.7 


4.03713x10" 


6 


5.23450x10" 


-6 


2 


0.666667 


-0 


0888889 


0.8 


4.20408x10" 


6 


5.77575x10" 


-6 


2 


0.666667 


-0 


0888889 


0.9 


4.59416x10" 


6 


6.26060x10" 


-6 


2 


0.666667 


-0 


0888889 



Table3: Errors and invariants for 
n = 100, a = -15, b = 15, St = 0.001, T = 0.1, . . . , 0.9. 



Time 


Loo 


L 2 


h 


h 


h 


0.01 


7.90227x10" 


-V 


4.72378x10" 


-V 


2 


0.666667 


-0 


0888889 


0.02 


1.03034x10" 


-6 


7.25023x10" 


-7 


2 


0.666667 


-0 


0888889 


0.03 


1.14705x10" 


-6 


9.14997x10" 


-7 


2 


0.666667 


-0 


0888889 


0.04 


1.21740x10" 


-6 


1.07375x10" 


-6 


2 


0.666667 


-0 


0888889 


0.05 


1.28565x10" 


-6 


1.21327x10" 


-6 


2 


0.666667 


-0 


0888889 


0.06 


1.41361x10" 


-6 


1.33912x10" 


-6 


2 


0.666667 


-0 


0888889 


0.07 


1.51656x10" 


-6 


1.45490x10" 


-6 


2 


0.666667 


-0 


0888889 


0.08 


1.60373x10" 


-6 


1.56261x10" 


-6 


2 


0.666667 


-0 


0888889 


0.09 


1.67755x10" 


-6 


1.66392x10" 


-6 


2 


0.666667 


-0 


0888889 



Tablc4: Errors and invariants for 
n = 100, a = -15, b = 15, St = 0.001, T = 0.1, . . . , 0.9andT = 0.001, . . . , 0.009. 
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6.2 KdVB equation 



By considering /i = 0.1, e = 2, v = 0.005, 9 = 0.5, in Table 5, two kinds of error 
is calculated for n = 100, a = -100, b = 100, St = 0.02, T = 1, . . . , 9. 

Table 6 indicates errors for [i = 0.001, e = 1, u = 0.001, 9 = 0.5 and n = 
100, a = -40, b = 100, St = 0.05, T = 1, . . . , 9. 

In Table 7, error is obtained for n = 16, a = 0, b = 100, St = 0.00001, T = 
0.0001, . . . , 0.0009 and fi = 0.001, e = l,v= 0.001, 6 = 0.5. 

Table 8 shows errors for n = 16, a = -40, 6 = 40, St = 0.02, T = 1, . . . , 9 
and ^ = 0.1, e = 2, v = 0.005, (9 = 0.5. 

In Table 9, errors are calculated for n — 16, a — 8, b = 99, St = 0.05, T = 
1, . . . , 9 and [i — 0.001, e = 1, v = 0.001, 9 = 0.5. Errors reduces by decreasing 
of time steps. 

Figure 3, indicates solutions of KdVB equation for n = 100, a = —100, b = 
100, St = 0.02, T = 1 and fi = 0.1, e = 2,v = 0.005, 9 = 0.5. 

Figure 4, indicates solutions of KdVB equation for n = 100, a = —40, b = 
100, St = 0.05, T = 1 and fi = 0.1, e = l,v = 0.1, (9 = 0.5. 

Figure 5 shows solutions of KdVB equation for n = 100, a = —40, b = 
100, St = 0.05, T = 1 and fi = 0.01, e = 1, v = 0.01, 9 = 0.5. 

In Figure 6, numerical solution of KdVB equation is shown for n — 100, a — 
-40, b = 100, St = 0.05, T = 1 and fi = 0.001, £ = 1, v = 0.001, (9 = 0.5 



Time 




L 2 


1 


2.06127x10" 


-0 


4.62479x10" 


-b 


2 


1.03034x10" 


-6 


7.25023x10" 


-7 


3 


1.14705x10" 


-6 


9.14997x10" 


-7 


4 


1.21740x10" 


-6 


1.07375x10" 


-6 


5 


1.28565x10" 


-6 


1.21327x10" 


-6 


6 


1.41361x10" 


-6 


1.33912x10" 


-6 


7 


1.51656x10" 


-6 


1.45490x10" 


-6 


8 


1.60373x10" 


-6 


1.56261x10" 


-6 


9 


1.67755x10" 


-6 


1.66392x10" 


-6 



Table5: Errors for /i = 0.1, e = 2, v = 0.005, 9 = 0.5 and 
n = 100, a = -100, b = 100, St = 0.02, T = 1, . . . , 9. 



Time 




Li 


1 


5.49372x10" 


-7 


1.00098x10" 


-b 


2 


1.09951x10" 


-6 


1.99828x10" 


-6 


3 


1.65039x10" 


-6 


2.99190x10" 


-6 


4 


2.20200x10" 


-6 


3.98187x10" 


-6 


5 


2.75434x10" 


-6 


4.96819x10" 


-6 


6 


3.30738x10" 


-6 


5.95089x10" 


-6 


7 


3.86111x10" 


-6 


6.92997x10" 


-6 


8 


4.41551x10" 


-6 


7.90547x10" 


-6 


9 


4.97059x10" 


-6 


8.87738x10" 


-6 
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Table6: Errors for fj, = 0.001, e=l,v = 0.001, 6 = 0.5 and 
n = 100, a = -40, b = 100, St = 0.05, T = 1, . . . , 9. 



Time 


-'-'oo 


L 2 


1 


1.84693x10" 


12 


6.66561x10" 


12 


2 


3.69387x10" 


12 


1.33287x10" 


11 


3 


5.54080x10" 


12 


1.99924x10" 


11 


4 


7.38774x10" 


12 


2.66562x10" 


11 


5 


9.23467x10" 


12 


3.33201x10" 


11 


6 


1.10816x10" 


11 


3.99840x10" 


11 


7 


1.29285x10" 


11 


4.66479x10" 


11 


8 


1.47755x10" 


11 


5.33118x10" 


11 


9 


1.66224x10" 


11 


5.99758x10" 


11 



Table7: Errors for fj, = 0.001, e=l,v = 0.001, 9 = 0.5 and 
n = 16, a = 0, b = 100, St = 0.00001, T = 0.0001, . . . , 0.0009. 



Time 




£ 2 


1 


7.82718x10" 


-8 


3.31435x10" 


-V 


2 


1.57051x10" 


-7 


6.62416x10" 


-7 


3 


2.36335x10" 


-7 


9.92942x10" 


-7 


4 


3.16123x10" 


-7 


1.32301x10" 


-6 


5 


3.96412x10" 


-7 


1.65261x10" 


-6 


6 


4.77201x10" 


-7 


1.98176x10" 


-6 


7 


5.58486x10" 


-7 


2.31043x10" 


-6 


8 


6.40266x10" 


-7 


2.63863x10" 


-6 


9 


7.22539x10" 




2.96637x10" 


-6 



Table8: Errors for /i = 0.1, e = 2, v = 0.005, = 0.5 and 
n = 16, a = -40, 6 = 40, St = 0.02, T = 1, . . . , 9. 



Time 


Loo 


£ 2 


1 


2.26009x10" 


-8 


8.83534x10" 


-8 


2 


4.52028x10" 


-8 


1.76687x10" 


-7 


3 


6.78055x10" 


-8 


2.65005x10" 


-7 


4 


9.04092x10" 


-8 


3.53306x10" 


-7 


5 


1.13014x10" 


-7 


4.41592x10" 


-7 


6 


1.35619x10" 


-7 


5.29861x10" 


-7 


7 


1.58226x10" 


-7 


6.18113x10" 


-7 


8 


1.80833x10" 


-7 


7.06350x10" 


-7 


9 


2.03441x10" 


-7 


7.94570x10" 


-7 



Table9: Errors for \i = 0.001, e=l,v = 0.001, = 0.5 and 
n = 16, a = 8, 6 = 99, £t = 0.05, T = 1, . . . , 9. 
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u 




Fig3: Solutions of KdVB equation for n = 100, a = -100, b 
0.02, T = 1 and fi = 0.1, s = 2, v = 0.005, 9 = 0.5 



u 




Fig4: solutions of KdVB equation for n = 100, a = —40, 6 
0.05, T = 1 and fj, = 0.1, e = 1, v = 0.1, 6» = 0.5. 



14 



u 




-0.005 L 

Fig5:solutions of KdVB equation for n = 100, a = -40, b = 100, St = 
0.05, T = 1 and fi = 0.01, s=l,v = 0.01, 9 = 0.5. 




Fig6: Solution of KdVB equation for n = 100, a = -40, b = 100, St = 
0.05, T = 1 and M = 0.001, £ = 1, v = 0.001, = 0.5 
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Conclusion 



The collocation method using sine basis is applied for solving KdV and KdVB's 
equations. Three invariant of conservation laws are calculated for KdV equation. 
The method is computationally attractive and results are shown through tables 
and figures. 
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